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ABSTRACT 

We measure black hole masses for 71 BL Lac objects from the Sloan Digital Sky Survey 
(SDSS) with redshifts out to z ^ 0.4. We perform spectral decompositions of their 
nuclei from their host galaxies and measure their stellar velocity dispersions. Black 
hole masses are then derived from the black hole mass - stellar velocity dispersion 
relation. We find BL Lac objects host black holes of similar masses, ~10*'^ Mq, with 
a dispersion of ±0.4 dcx, similar to the uncertainties on each black hole measurement. 
Therefore, all BL Lac objects in our sample have the same indistinguishable black 
hole mass. These 71 BL Lac objects follow the black hole mass - bulge luminosity 
relation, and their narrow range of host galaxy luminosities confirm previous claims 
that BL Lac host galaxies can be treated as standard candles. We conclude that the 
observed diversity in the shapes of BL Lac object spectral energy distributions is not 
strongly driven by black hole mass or host galaxy properties. 

Key words: galaxies: active - BL Lacertae objects: general - galaxies: jets - galaxies: 
kinematics and dynamics 
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1 INTRODUCTION 

BL Lac objects are a very rare class of Active Galactic 
Nuclei (AGNs) thought to be unified with Fanaroff-Riley I 
(FR I, iFanaroff fc Rilevlll974l 'l radio gala xies viewed nearly 
along the axis of a relativistic jet (e.g., see lBlandford fc ReesI 
ll978l : IUrrv fc Padovanilll995l '). BL Lac objects are among the 
brightest radio. X-ray, and gamma-ray sources in the entire 
sky, they are highly polarized (>2-3%), and their multiwave- 
length emission is variabl e on timescales ranging from min- 
utes to decades (e.g., see lKollgaardlll9S)3 : lUrrv fc Padovanil 
[l99i). Perhaps the most striking property of BL Lac objects 
is their optical spectra, which are devoid of strong emission 
lines because radiation from the Doppler boosted jet is so 
dominant. 

Black hole mass measurements for a large number of 
BL Lac objects would have a wide range of applications 
- from investigations aimed toward a better understanding 
of BL Lac statistical properties, to studies focusing more 
generally on relativistic jets. For example, AGNs are pow- 
ered by accretion onto a supermassive black hole, making 
black hole mass a fundamental parameter - i.e., it sets an 
AGN's upper luminosity limit, and it determines the tem- 
perature of the accretion disc. Also, BL Lac spectral en- 
ergy distributions (SEDs) are dominated by jet emission, so 
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they are excellent probes of relativistic jet physics. BL Lac 
objects show impressively diverse SEDs from object to ob- 
ject, with a synchrotron cutoff occurring anywhere from 
the near-infrared to the soft-X-ray wavebands (e.g., see 
IPadovani fc Giommil Il995h . A large sample of black hole 
masses can thus be used to investigate if black hole mass 
or anoth er parameter (s) is the primar y driver of a jet's SED 
(e.g., see lGhisellini fc TavecchiclbOOSi '). Even more generally, 
a sample of BL Lac objects with well determined black hole 
masses could be incorporated in studies on the scaling of 
relativistic jets with mass, through, for example, compari- 
son with X-ray binaries in the hard sp ectral state (e.g., see 
iMerloni et aLH2003l : iFalcke et aLll2004l l . 

The weak-lined nature of BL Lac spectra, however, in- 
troduces unique challenges for obtaining reliable black hole 
mass measurements. For other classes of AGN, black hole 
masses can be derived by measuring the widths of gas emis- 
sion lines and by measuring the continuum luminosity - 
here, one assumes the gas in the broad emission line re- 
gion is virialized and adopts a relation between continuum 
lu minosity and size o f the broad emission line region (e.g., 
see iKaspi et al.ll2000l ). However, this technique cannot be 
employed for AGN with featureless spectra. Even if weak 
emission is detected (and if one assumes correlations be- 
tween omission line strength and black hole mass extend to 
weak emission) , the resultant BL Lac black hole masses will 
be unreliable unless the continuum luminosity is corrected 
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for Doppler boosting. Other techniques for weighing BL Lac 
central black holes should be seeked. 

An alternative method utilizes the dynamics of the 
stars near a galaxy's center, i.e., 'dynamical' mass estimates. 
This is done by measuring the line-of-sight stellar velocity 
dispersion (a) from the widths of stellar absorption fea- 
tures, which are Doppler broadened due to the stars' mo- 
tions. Thus, if one can decompose an AGN's host galaxy 
from the central point source and measure a, then black 
hole masses can be estimated f rom empirical correlations 
(the Mbh-o' re l ation , e.g., see iFerrarese &: MerrittI I2OO0I : 



iGebhardt et al.l I2OO0I : iTremaine et all I2OO2I '). Black hole 
masses can also be estimated from the host galaxy's lumi- 
nosit y (e.g.. iKormendv fc Richstone|[l995l : iMagorrian et al.l 
Il998l ). and sizable samples of BL Lac black hole masses 
(^--^60- 120 objects) have indeed been derived in such a way 
(e.g., IWu et all [2OO9I : IXu et al.l [2OO9I '). However, the black 
hole mass - galaxy luminosity correlation is not as tight as 
the Mbh-o" relation, and black hole masses derived from 
galaxy luminosities can sometimes be overestimated by up 
to tw o orders of magnitude (e.g., see IFerrarese fc MerrittI 

bood) . 

Given the above considerations, we prefer to estimate 
black hole masses using stellar velocity dispersions. How- 
ever, in the past, there was not an abundant supply of 
BL Lac objects with high-quality spectra with enough ap- 
parent host galaxy light to perform the requisite measure- 
ments of a. For this reason, subsets of BL Lac objects 
with dynamical black hole estimates tend to contain at 
most ~30 objects (e.g., see iFalomo et al" I I2OO3I : IWoo et al.l 
I2OO5I ). Fortunately, relatively large BL Lac samples con- 
taining hundreds of objects have recently come online (e.g., 
iTurriziani et al.ll2007l : lMassaro et al.ll2009l ). Here, we use the 
723 object BL Lac sam ple assembled fro m t he Sloan Digi- 
tal Sk y Survey (SDSS. lYork et al.ll2000l ) bv iPlotkin et al.l 
l|2010l . hereafter PIO). Each object has a high-quality SDSS 
spectrum, and the sheer size of this sample ensures a rela- 
tively large number of BL Lac spectra show a strong enough 
host galaxy component to meas ure a . We attempt black 
hole mass measurements for 143 IPIOI BL Lac objects, for 
which we derive 71 black hole masses. iLeon-Tavares et al.l 
l|2010l ) recently presented a similarly sized sample of dynam- 
ical blac k hole masses from our older radio-selected BL Lac 
catalog (|Plotkin et al.l [20081 ). In iJHwe describe the parent 
BL Lac sample, and we explain our algorithm for estimat- 
ing black hole masses in ^ The measured distribution of 
black hole masses and host galaxy luminosities are discussed 
in ij4l and we comment on the import of black hole mass 
in determining the shape of BL Lac SEDs. Finally, our 
conclusions are summarized in SJS] Throughout, we adopt 
Ho = 71 km s"^ Mpc"^ SI™ = 0.27, and = 0.73. 



2 THE BL LAC SAMPLE 



The IPIOI optically selected BL Lac sample contains 723 ob- 
jects from 8250 deg^ of SDSS spectroscopy. About 75% of 
these BL Lac objects were discovered by the SDSS. In or- 
der for a source to be retained as a BL Lac candidate, its 
SDSS spectrum cannot show any emission line with rest- 
frame equivalent width (REW) stronger than 5 A, and its 
Ca II H/K depression must be smaller than 40%. There is 



no explicit requirement for an object to be a radio or X-ray 
emitter, but the sample was correlate d post-selection with 
the NRAO VLA Sky Survey (NVSS. [Coiidon et al.|[l998l ) 
and with the Faint Images of the Radio Sky at Twenty-cm 
(FIRST. [Becker et al.|[l995l ) radio surveys, as well as with 
the ROSAT All Sky Survey (RASS lVoges et al.lll999l . I2OO0I ) 
in the X-ray. We refer the reader to Fig for details. 

About one-third of SDSS BL Lac objects have spectro- 
scopic redshifts derived from host galaxy spectral features 
in their SDSS spectra. Here, we limit ourselves t o 143 spec- 
tra that have reliable spectroscopic redshifts in jPlOl . that 
match to a FIRST and/or NVSS radio source and are radio- 
loud (i.e., radio-to-optical flux ratio is lar ger than 10, see 



iKellermann et aL [ 1 19891 : IStocke et al.1ll992l )Fl and that have 
z < 0.4. Objects with z > 0.4 are too distant to show sig- 
nificant host galaxy radiation in their SDSS spectra. We 
attempt to measure black hole masses for each of these 
143 BL Lac objects. These 143 objects are among the most 
weakly beamed BL Lac objects in [PlOl (or else they would 
not show a host galaxy component in their optical spectra) , 
thereby minimizing orientation effects. 



3 BLACK HOLE MASS MEASUREMENTS 

To measure black hole masses we first remove the contribu- 
tion from the AGN to the observed spectrum ( i]3.H) . and we 
then measure the line of sight stellar velocity dispersion (a) 
in each decomposed host galaxy spectrum ( >i3.2|l . We assess 
the reliability of the spectral decompositions in ^3.31 and in 
that section we identify 72 spectra with unreliable decom- 
positions (most have poor fits based on their reduced x^)- 
Finally, we estimate black hole masses and host galaxy lumi- 
nosities for the other 71 spectra ( H3.4l and l3.5l respectively). 



3.1 Spectral Decomposition 

We assume each SDSS spectrum is a combination of two 
components: the thermal contribution from the host galaxy, 
and the non-thermal emission from the relativistic jet. We 
model the latter as a power law: = fi,g {v / vo)'"" , where 
a^j is the spectral index and f^o is the flux density at refer- 
ence frequency vq. We choose !/o=c/(6165 A), which is the 
reference frequency of the SDSS r fllter. To find the best fit 
power-law and host galaxy parameters, we make slight mod- 
ifications to the Gas And Absorption Line Fitting algorithm 
(GANDALF)[3 as described below. 

GANDALF iteratively fits an observed galaxy spectrum 
with combinations of stellar templates convolved by the 
best-fit line of sight velocity dispersion; it also simultane- 
ously fits a user defined list of gas emission lines, model- 
ing the lines with Gaussian templates. The Penalized Pixel- 



^ IPIOI include a population of ~10'^ weak-lined AGN lacking 
strong radio emission, many of which are unlikely best unified 
with BL Lac objects. This constraint on radio emission therefore 
ensures wc consider only normal BL Lac objects. 

GANDALF was developed by the SAURON team and is avail- 
able from http : //w ww. strw. leldenuniv .nl/sauron| Also see 
ISarzi et al.i (,2006.) . 
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Fitting method (pPXF0, see ^3.2[l of lCappellari fc Emsellenil 
(|2004l ) is used by GANDALF to provide initial guesses for the 
stellar continuum fits and velocity dispersions. Nonlinear 
fits are then performed with the Bounded- Variables Least- 
Squares algorithm (BVL^j)) a^nd with Levenberg-M arquardt 
least-squares minimization (see iMarkwardtl |2009| ) . An im- 
portant feature of GANDALF is that it fits the stellar contin- 
uum and gas emission simultaneously, instead of masking 
spectral regions containing emission lines from the contin- 
uum fits. Thi s technique is like ly to yield less biased stellar 
continua (see lSarzi et ahllioOfih . 

We modified GANDALF to also include a power law com- 
ponent from the AGN, and we then ran the modified code on 
the 143 BL Lac objects showing a host galaxy component 
in their SDSS spectra. All fits are performed in the rest- 
frame. The power law is constrained to always be positive 
and to not account for more than 80% of the observed flux 
at 3900 A. The latter constraint is to ensure the power law 
is physically meaningful: if the power law accounted for all 
of the observed flux, then we would not see any contribution 
from the host galaxy to the BL Lac spectrum . We use stel- 
lar population models jVazdekis et al.ll2010l ) based on the 
Medium-resolution Isaa c Newton Telescope library o f em- 
pirical spectra (MILES. ISanchez-Blazguez et al.]|2006l ). and 
the SDSS BL Lac spectra are convolved to match the spec- 
tral resolution of the MILES templates (81 km s^^). Because 
BL Lac emission lines are by deflnition very weak, we only 
flt gas emission for [O II], [O III], H/3, Ha, [N II], and 
[S II], and we exclude spectral regions containing Oxygen 
sky lines from the flts. Sample spectral decompositions are 
shown in Fig. [T] 



3.2 Velocity Dispersion Measurements 

For each BL Lac spectrum, we subtract the best flt power 
law and emission lines, yielding a spectrum of just the host 
galaxy. We then measure the host galaxy's velocity disper- 
sion, a, with the pPXF method. The pPXF algorithm flts 
galaxy spectra by convolving sets of stellar templates with a 
best-fit velocity dispersion. We again use the MILES stellar 
templates, and we adopt the best-fit combination of stellar 
templates found during the GANDALF spectral decomposition 
in t]3.1l Since we remove any (weak) emission lines from our 
host galaxy spectra, there is no need to mask out regions of 
the spectra containing gas emission. We do, however, still 
omit regions covering Oxygen sky lines from the pPXF fits. 

We estimate errors on a via Monte Carlo simulations 
(see ICappellari fc Emsellemll2004h . For each pixel's flux den- 
sity, /A,i, we add noise ±ei, where is drawn randomly 
from a normal distribution with standard deviation equal 
to each pixel's flux density measurement error. We then re- 
measure the velocity dispersion with pPXF, repeating 100 
times. We adopt the standard deviation in the 100 veloc- 
ity dispersion measurements as the formal uncertainty on a, 
which is typically better than 10%. Uncertainties on power 



IDL code for running pPXF can be down- 
loaded from Michele Cappcllari's web page: 
|http : //www- astr o .physics ■ ox ■ ac.uk/-mxc/idl/ 
* Code is also available from Michele Cappcllari's web page: 
|http : //www- astro .physics . ox . ac.uk/~mxc/idl/ 
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Figure 1. Three sample SDSS spectra, in their rest-frames, con- 
volved to the MILES spectral resolution, shown as black solid 
lines. The best GANDALF fits are overdrawn as red solid lines, and 
the red dashed lines illustrate the power law components. This 
figure appears in color in the online version of this article. 



law indices are similarly derived by randomly adding noise 
and re-running GANDALF 100 times per spectrum. The best- 
flt stellar velocity dispersions and their uncertainties are in- 
cluded in Table |T] 

In Table[T]we also include the best flt power law indices, 
a^, and the fraction of the observed flux density accounted 
for by the power law at rest-frame 6165 A, fpi. The measured 
values of are typical for BL Lac objects, with an average 
(qi,) = 0.70 ± 0.56. The recovered values of fpi (with an 
average value 0.34 ± 0.12) are as expected. If the power law 
is too weak it would not be a BL Lac object, and if the 
power law is too strong the host galaxy would not be seen 
in the optical spectrum. The fraction of light accounted for 
by the power law is systematically higher if the calculation 
is instead done at a wavelength bluer than 6165 A. 



3.3 Assessing the Quality of the Decomposition 

There are several reasons why some of the spectral decompo- 
sitions might fail. All 143 SDSS BL Lac spectra are of high 
enough quality to derive redshifts from host galaxy spec- 
tral features, and to allow reliable BL Lac classiflcations. 
However, the amount of the host galaxy flux for the most 
highly beamed objects in this subset is not strong enough 
to obtain accurate spectral decompositions, and some other 
objects do not have high enough signal-to-noise. We thus ex- 
clude 54 objects with reduced Xr > 2 (each GANDLAF flt has 
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Table 1. Black Hole Masses and Other Measured Parameters 
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~4200 degrees of freedom) , and we remove an additional 12 
objects for which the spectral decompositions did not con- 
verge, or the error on the velocity dispersion is larger than 
40 km s~^ (i.e., half the spectral resolution of the MILES 
stellar templates). Upon visual examination of each spec- 
tral fit, we remove an additional six objects with poor fits. 
Most of these six objects have Ca II H/K break strengths 
smaller than 15%, and their continua are too dominated by 
the power law to obtain reliable a measurements. We are left 
with 71 spectra for which we determine black hole masses. 
For the remainder of this paper we only consider these 71 
objects. 

The 71 successfully decomposed objects are clearly bi- 
ased toward those with strong host galaxy components. The 
amount of stellar flux in BL Lac spectra depends primarily 
on orientation. That is, mor e weakly beamed BL Lac objects 
show more stellar light (e.g.. lLandt et al" ]|2002l : |Plotkin" et al.l 
l2008h . and our 71 object subset is therefore biased toward 
larger viewing angles. All of the 71 BL Lac objects with 
successful decompositions have Ca II H/K breaks stronger 
than 15%, roughly translating to viewing angles larger than 
about 15-25 deg (see Fig. 8 of lLandt et"al]|2002l ). Since the 
above bias is primarily geometric, we make the reasonable 
assumption that our analysis and conclusions on this subset 
can be extrapolated to the entire BL La c population. 

As pointed out by I Woo et al.l |200^, less massive black 
holes live in less luminous galaxies and are therefore more 
difficult to measure since the host galaxy is harder to see, 
especially at increasing redshift. This selection effect could 
artificially narrow the observed stellar velocity (and there- 
fore black hole mass) distribution. We however do not see 
any correlation between velocity dispersion and redshift in 
our sample of 71 objects, and there is no obvious deficit of 
low mass black holes at high-redshift. We therefore do not 
believe this strongly affects our sample. Regardless, since we 
do not attempt to derive a complete sample of BL Lac ob- 
jects with black hole masses here, this potential bias will not 
alter our conclusions. 



3.4 Measuring Black Hole Masses 

The Mbh-o" relation is parameterized as 

log(MBH/M0) =aH-/31og(o-/cro), (1) 

with o-Q = 200 km s~ \ We use the coefficients derived by 
iTremaine et all (|2002l ) from a sample of 31 nearby galaxies: 
a = 8.13 ± 0.06 and /3 = 4.02 ± 0.32. To estimate error bars, 
we propagate the uncertainties on our velocity dispersion 
measurements and the errors on a and p. We add that un- 
certainty in quadrature with the Mbh-o" re lation's intrinsic 
scatter in logMen: o"int ~ ±0.30 dex fsee ITremaine et al.l 
|2002| ). Typical error bars on our black hole measurements, 
including the intrinsic scatter of the Mbh — o" relation, are 
±0.30-0.35 dex. The measured black hole mass distribution 
is shown in Fig[2l and masses are listed along with the other 
fitted parameters (i.e., power law indices, velocity disper- 
sions, etc.) in Table [T] 

There appears to be a small number of black hole mass 
measurements that deviate from the mean black hole mass 
in Figure [2] (i.e., log Mbh < 7.5 and logAfaH > 9.0). We 
re-examined their spectral decompositions visually, and we 
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Figure 2. Distribution of black hole mass measurements, derived 
with Equation [T] 



find no compelling reason to remove them from the sample. 
Instead, we attribute those points to populating the high 
and low sigma tails of the mass distribution, and to small 
number statistics. A Kolmogorov-Smirnoff test shows our 
measured black hole mass distribution is consistent with a 
Gaussian, in which case we would expect 3.6 black hole mass 
measurements to deviate from the mean by ±2cr. Indeed, 
four black hole mass measurements populate the 2a tail. 



3.5 Measuring Host Galaxy Luminosities 

After performing the spectral decomposition we are able to 
construct spectra of each BL Lac object's host galaxy. We 
use these spectra to estimate the absolute magnitude each 
host galaxy would have in the Cousins Rc filter at z = 0. 
First, we measure the host galaxy fiux at rest-frame 6581 A 
(the effective wavelength of the Rc filter) , and we then apply 
a bandpass correction and correct for extinction from the 
Milky Way (using the dust maps of lSchlegel et al.|[l998l ). 

The 3" apertures of each SDSS spectroscopic fiber does 
not include all of the host galaxy's fiux, and an aperture 
correction must be applied. To determine the fraction of 
light covered by each spectroscopic fiber, we perform an im- 
age decomposition on ea ch source's sdss i filter image us- 
ing the GALFIT software l|Peng et al.ll2002l ). We model each 
BL Lac object with a point source^ and a de Vaucouleur 
profile (BL Lac objects prob ably live exclusively in ellipti- 
cal galaxies. [Urrv et aLllioOol ). We are able to measure the 
half light effective radius re for 29 BL Lac hosts, requiring 
re > 3 pix, and that the best-fit host galaxy and central 
point source models are centered on the same coordinates. 
We adopt these values to calculate the aperture correction 
for these 29 objects. For the remaining objects, we assume 
re — 10 kpc, the median effective radius for the 29 objects. 



^ The point-spread function at the location of each 
BL Lac object on the SDSS camera is extracted with 
the read_PSF code provided by the SDSS collaboration: 
http: / /www. sdss . org/dr7/products/ images/read_psf .html 
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This value is also consistent with the median effective ra- 
dius from other BL Lac host galaxy studies (e.g-. IUrrv et alj 
I2OOOI ). 

The above aperture correction contributes the largest 
source of uncertainty to our luminosity estimates. For the 29 
objects with re measurements, the uncertainties on re are 
typically of the order ±1 kpc (0.6 pix), introducing typical 
errors on the host galaxy absolute magnitudes of ±0.06 mag. 
The systematic error introduced for the other 42 objects by 
assuming re = 10 kpc is more severe. To estimate the un- 
certainty, we calculate the host galaxy luminosity for the 
29 objects with measurements by assuming = 10 kpc, 
and we compare to the host galaxy luminosity estimates 
using each re measurement. There is no systematic offset 
between the two absolute magnitude measurements; how- 
ever, the measurements disagree with a standard deviation 
of ±0.3 mag. We therefore adopt a systematic uncertainty 
of ±0.3 mag for the 42 objects without re measurements. 
Note, the error's dependence on redshift is weak for z > 0.1. 

We add the above uncertainties in quadrature to the 
statistical errors estimated from Monte Carlo simulations, 
where we randomly added noise to each spectrum and re- 
ran each spectral decomposition 100 times (as described in 
i)3.2p . These uncertainties, including error from the aperture 
corrections, are listed in Table [1] For the 29 objects with 
re measurements, we do not estimate absolute magnitudes 
using the total integrated host galaxy apparent magnitude 
from their best-fit de Vaucouleur profiles because the re- 
quired K-corrections would introduce additional uncertain- 
ties. 




log (L/L,) 

Figure 3. Logarithm of black hole mass vs. logarithm of the ratio 
of radio (5 GHz) to X-ray (1 keV) specific luminosities. The lat- 
ter parameterizes the shape of each object's SED. Error bars are 
omitted for clarity - a typical (Icr) black hole mass measurement 
error is shown in the top left corner. The measurement errors 
on log (Lr/Z/x) are smaller than the data points. We do not ob- 
serve a significant correlation indicating one or more parameters 
other than black hole mass are important for determining the 
SED shape. 



4 DISCUSSION 

4.1 Black Hole Masses 

We measure an average black hole mass of (log Mbh) = 
8.54 ± 0.047 Mq (the quoted uncertainty is the error of the 
mean), with a standard deviation of ±0.40 dex. The disper- 
sion in the distribution of black hole masses is comparable to 
the typical measurement error of ±0.30-0.35 dex. Thus, it is 
possible that the observed spread of black hole masses can 
be attributed only to measurement error, and we conclude 
BL Lac black hole masses are virtually indistinguishable. 

Our black hole estimates are consistent with those 
in the literature. For example, iFalomo et all (|2003l l find 
(log Mbh) = 8.57 Mq for a sample 1 2 BL Lac objects with 
dynamical black hole mass measures. IWoo et al.1 (|2005l ) find 
32 BL Lac objects with dynamical black hole mass mea- 
sures to have a dispersion from 4 x lO'^ — 6 x 10* ^©i S'lid 
(log Mbh) ~ 8.30 Mq. They note, however, that their sam- 
ple is relatively shallow {{z) ~ 0.17), and it thus may not 
probe a large enough volume to recover a substantial num- 
ber of massive central black holes. Indeed, when they in- 
clude an additional ~30 BL Lac objects with black hole 
mass measures based on host galaxy luminosity, their ~60 
object sample has (z) ^ 0.31 and a dispersion in black hole 
mass measures that extends to 4 x 10^ Mq (identical to the 
largest black hole mass measured in our 71 object sample). 

It is interesting that the black hole mass distribution is 
so narrow, because BL Lac SEDs can be incredibly diverse: 
BL Lac objects are often classified as high-energy peaked 



(HBL) or low-energy peaked (LBL) BL Lac objects, depend- 
ing on the cutoff frequency o f their synchrotr o n rad iation 
fsee lPadovani fc Giomrnilll995l 'l. lNieppola et al.1 (j2006l 'l show 
BL Lac synchrotron radiation cutoff frequencies can span 
close to ten orders of magnitude from object to object. Yet, 
we find that differences in SED shapes must not significantly 
depend on black hole mass. 

We parameterize the SED shape of our 71 BL Lac 
objects with \og{Lr / Lx), where Lr and Lx are the spe- 
cific luminosities (in erg s~^ Hz"^) at rest-frames 5 GHz 
and 1 keV respectively, and their values are taken from 
iPlOl . The SEDs of BL Lac objects with smaller Lr/Lx 
values ha ve synchrotro n cuto ffs at higher frequencies. 
We follow IPlotkin et al.l ()2008h and classify objects with 
log (L^/Lx) < 5.76 as HBLs. Our sample of 71 BL Lac 
objects includes 59 HBLs and 12 LBLs; the relatively larger 
number of H BLs in our sample is an artifact of SDSS BL Lac 
selection (see lPlol ). A Kolmogorov-Smirnoff test indicates a 
p=0.47 chance that the HBL and LBL black hole mass dis- 
tributions come from the same population, and their black 
hole mass distributions are therefore not statistically differ- 
ent. 

The logarithm of black hole mass vs. \og{Lr / Lx) is 
shown in Fig. [3] and we see no statistically significant cor- 
relation. The linear Pearson correlation coefficient between 
log(Lr/Li,) and log Mbh is r = 0.018, with a probability 
p — 0.44 of randomly finding a stronger correlation, treat- 
ing Lr/Lx limits as exact. We find r — 0.036 and p — 0.40 
if we only consider the 55 objects with X-ray detections in 
RASS. Thus, differences between BL Lac SEDs are likely 
not driven solely by black hole mass. This is consistent with 
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iGhisellini fc Tavecchicl (|2008l ). who show the observed ranee 
of BL Lac synchrotron cutoff frequencies can be replicated 
by assuming all BL Lac objects have the same black hole 
masses and varying only their Eddington normalized accre- 
tion rates (see their Fig. 5). We do not derive normalized ac- 
cretion rates for the 71 BL Lac objects with black hole mass 
measurements here because, without knowing their Doppler 
factors, we cannot accurately debeam their luminosities. 

4.2 Host Galaxies 

A correlation between black hole mass and host galaxy ab- 
solute magnitude is found (Fig. [4|, confirming that 
BL Lac object s follow the black hole mass - bulge rela- 
tionship (e.g., iKormendv fc Richstond Il995l ). Our sample 
contains the largest number of BL Lac objects for which 
the correlation has been tested . We p erform a linear re- 
gression using the technique of iKellvl (|2007l ') , which takes 
a Bayesian approach and accurately accounts for measure- 
ment errors in both variables. We find log [Mbh /Mq^ = 
(-4.12 ± 2.99) - (0.55 ± 0.13) Mr^, and a linear Pearson's 
correlation coefficient r = —0.513 (p < 10~^ chance of ran- 
domly finding a stronger corr elation). Our regress i on ag rees 
with the correlation found bv lMcLure fc Dunlod (|200a ) for 
a sample of 90 active and inactive galaxies covering similar 
redshift ranges and host galaxy luminosities (dashed line in 
Fig.©. 

There are a small number of points that deviate from 
the black hole mass - bulge relation in Fig. U (i.e., those with 
logMsH < 7.5 and logMsH > 9.0). There is no motivation 
to remove those objects from our sample based on the qual- 
ity of their spectral decompositions, and we thus include all 
objects in the above regression. We however note they merit 
further study. 

There is no reason to assume the observed correlation 
is driven by correlated errors between our BL Lac black 
hole mass and host galaxy luminosity measurements. We 
test this by randomly adding noise to each observed spec- 
trum and rerunning our spectral decomposition 100 times 
for each BL Lac object. No source shows a significant cor- 
relation between its 100 pairs of logMen and Mr^: the 
average Pearson's correlation coefficient for the 71 sources 
is (r) = —0.01, and no source shows an anti-correlation 
stronger than r = —0.2. 

That BL Lac objects obey the same black hole mass 
- bulge relation as other galaxies indicates that the black 
hole formation histories of BL Lac objects are similar to 
other AGN. This is further support that parameters other 
than black hole mass (such as orientation, accretion rate, 
etc.), must be important for dictating if an AGN looks like a 
BL Lac phenomenon, and then for determining if the BL Lac 
object is an LBL or HBL. 

With the black hole mass - bulge relation in mind, we 
expect a small scatter in the BL Lac host galaxy luminos- 
ity distribution, since they have a relatively narrow range 
of black hole masses. This is indeed observed, as shown in 
Fig. [5] The BL Lac host galaxies are well fit by a Gaus- 
sian with (Mr^) = —22.9 mag and a standard deviation of 
just ±0.4 mag. This co nfirms the "standard candle" result of 
ISbarufatti et all (|2005h . who showed 64 BL Lac objects with 
their host galaxies imaged by the Hubble Space Telescope 
have host galaxy absolute magnitudes well-fit by a Gaus- 
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Figure 4. Logarithm of black hole mass vs. host galaxy abso- 
lute magnitude in the Cousins Rc filter. Error bars on absolute 
magnitude are shown for the 29 objects with effective radius mea- 
surements. For clarity, a typical error bar for the other 42 objects 
(where we assume re = 10 kpc) is shown in the top left corner, 
as is the typical la uncertainty for all 71 black hole mass mea- 
surements. The solid red line s hows our be st-fit linear regression 
using the Bayesian method of iKellvl ll2007t) . and the red dotted 
lines illustrate the ±3cr uncertainties. The bl ack hole mass - bulge 
relationship from iMcLure &: Dunlod ||2002| ) is overplottcd as a 
black dashed line for comparison. This figure appears in color in 
the online version of this article. 



sian with mean —22.8 mag and standard deviation ±0.5 mag 
(using the same cosmology employed herein). For the red- 
shifts considered here, the distribution is sufficiently narrow 
to estimate redshifts better than Az ~ ±0.07-0.08 (e.g ., 
see ISbarufatti et al. I l2005l . iMeisner fc Romanil [2OI0I . [PIO ). 
iMeisner fc Romanil (|2o'Tol ) and lPld apply the standard can- 
dle result to show that BL Lac objects lacking spectroscopic 
redshifts tend to be more distant on average, which can have 
important implications for understanding BL Lac evolution. 

We find no significant difference between the distribu- 
tion of host galaxy luminosities for HBLs and LBLs. The 
same host galaxy absolute magnitude {{Mr^) — —22.9 mag) 
is measured for both HBLs and LBLs, and the Student's t- 
test indicates a probability p = 0.977 that both distributions 
have the same mean assuming their variances are similar 
(p = 0.973 if we allow for different variances). The Student's 
t-test assumes the data follow a normal distribution, which 
is justified by our observations (see Fig. [5]). Even if we drop 
this assumption, the non-parametric Mann- Whitney U-Test 
still indicates the sample means are not significantly differ- 
ent (at the p = 0.379 level). Finally, a Kolmogorov-Smirnoff 
test gives a p = 0.820 chance the samples are drawn from 
the same parent population. We thus conclude that the ob- 
served diversity in BL Lac SEDs is not caused by differences 
in host galaxy properties. 
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Figure 5. Absolute magnitude of each decomposed host galaxy 
in the Rc filter. The distribution of absolute magnitudes is well 
described by a Gaussian with mean -22.9 mag and standard de- 
viation ±0.4 mag. This confirms earlier claims that the luminosi- 
ties of BL Lac host galaxies is narrow enough to treat them as 
standard candles. We do not see any differences in the host lu- 
minosities of LBLs or HBLs, indicating the host galaxy does not 
play a significant role in determining their SEDs. 

5 CONCLUSIONS 

AGN very likely influenced the formation of galaxies and 
heira rchical structure in the universe (e.g.. Idi Matteo et al.l 
obtaining reliable AGN black hole mass estimates is 
therefore an important endeavor. We presented dynamical 
black hole mass measurements for 71 BL Lac objects that 
show enough host galaxy contribution to their SDSS spectra 
to measure their stellar velocity dispersions. 

We find SDSS BL Lac objects have similar black hole 
masses, with the dispersion in their distribution comparable 
to the measurement uncertainty (~0. 30-0. 35 dex). We do 
not see any trend between black hole mass and the ratio of 
radio to X-ray luminosity, which serves as a proxy for SED 
shape. Our conclusion is that the diverse SED shapes exhib- 
ited by BL Lac objects is not controlled solely by black hole 
mass; another para meter, perhaps accretion rat e, must be 
important (e.g., see iGhisellini fc Tavecchioll2008l ). Dynami- 
cal black hole mass measures can only be obtained for the 
most weakly beamed BL Lac objects, or else their optical 
spectra do not show a strong enough host galaxy component. 
Therefore, the 71 BL Lac objects with black hole measures 
have similar orientation angles, and relativistic beaming is 
likely secondary in determining the shape of their SEDs. 
However, ori entation effects may not be completely negligi- 
ble (e.g., see lNieppola et al.|[200t ). 

We also recover a correlation between black hole mass 
and host galaxy luminosity that is consistent with the black 
hole mass - bulge lumin osity relation in the literature (e.g., 
iMcLure fc Dunlod[200^ ). This indicates the formation his- 
tories of BL Lac central black holes is similar to other AGN; 
we do not see any evidence that SED shape depends on 
host galaxy properties either. Our host galaxy decomposi- 



tion also supports claims that the distribution of BL Lac 
host galaxy luminosities is narrow enough to treat them as 
standard candles. One can thus derive "imaging redshifts" 
for BL Lac objects lacking spectroscopic redshifts if their 
host galaxies can be resolved from the nucleus, and redshift 
lower limits can be estimated when host galaxies are not 
detected. 

The black hole masses presented here expand the types 
of applications of SDSS BL Lac catalogs. For example, 
relativistic outflows are observed not just from supermas- 
sive black holes but also from their much smaller Galactic 
(~10 Mq) counterparts. The discovery of the fundamental 
plane of black hole accretion, a correlation between radio 
luminosity. X-ray luminosity, and black hole mass, demon- 
strates a remarkable connection between accreting black 
holes over eight orders of magnitude of black hole mass (e.g., 
iMerloni et aLll2003l : iFalcke et aLll2004l '). Part of our motiva- 
tion for measuring BL Lac black hole masses is to investigate 
how well BL Lac objects fit onto the fundamental plane; this 
will be discussed in a future paper. BL Lac objects are an 
interesting AGN subclass for fundamental plane studies be- 
cause their non-thermal emission is dominated by the jet. 
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